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In a previous publication local gauge invariant geometric variables were introduced 
to describe the physical Hilbert space of Yang-Mills theory. In these variables, the electric 
energy involves the inverse of an operator which can generically have zero modes, and thus 



its calculation is subtle. In the present work, we resolve these subtleties by considering a 
small deformation in the definition of these variables, which in the end is removed. The 
case of spherical configurations of the gauge invariant variables is treated in detail, as well 
as the inclusion of infinitely heavy point color sources, and the expression for the associated 
electric field is found explicitly. These spherical geometries are seen to correspond to the 
spatial components of instanton configurations. The related geometries corresponding to 
Wu-Yang monopoles and merons are also identified. 
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1. Introduction 

A new formulation of nonabelian gauge theory has recently been given (hereafter 
referred to as L), which is founded on a geometrical basis and seeks to define a setting in 
which gauge symmetry is implemented exactly and manifestly, even under approximations 
of the dynamics. Here, we will present further developments of this formulation as well as 
an example showing how it may allow one to obtain new insights into the reasons why the 
Yang-Mills field theory has a mass gap and produces a long range confining interaction 
between massive colored sources. We have introduced this geometric basis in the Hamilto- 
nian formulation because in the Lagrangian path integral formalism for a nonabelian gauge 
theory, the local gauge invariance does not really manifest itself as a symmetry, but rather 
as a redundancy in the path integral measure. Any approximations in that formalism are 
likely to introduce gauge artifacts precisely because local gauge invariance is not acting as 
a symmetry. In the strong coupling region of the theory, an alternative and perhaps better 
procedure might be one that avoids such gauge artifacts. This is possible in the canonical 
Hamiltonian formulation, because in "temporal" gauge, Aq = 0, there is a remaining local 
gauge invariance restricted to space-dependent transformations at a fixed time, and one 
can achieve the goal of treating this gauge invariance like a true quantum mechanical sym- 
metry. The local generators of such gauge transformations form an algebra and represent 
symmetries of the Hamiltonian which can be maintained exactly, even when the dynamics 
is done approximately. In view of this, one can understand better the statement that gauge 
invariance is not a symmetry away from a fixed-time formalism, and why approximations 
introduce gauge artifacts: the full classical gauge group includes time-dependent gauge 
transformations, and these are coupled to the dynamics. An approximate version of the 
dynamics is then likely to destroy any attempt to keep the full gauge invariance. 

We implement the gauge symmetry by considering a change of variables in the Hilbert 
space such that any function of the transformed variables is a singlet under the gauge group, 
i.e., is gauge invariant. In fact, this can be done in different ways [0-0. In I., we have 
introduced such a transformation of variables, and here we will develop this formalism in 
more detail. Our basic procedure is straightforward and simple to state: rather than using 
the space components of the vector potential, Af (a is a color index, i = 1,2,3 a space 
index), as fundamental coordinates in the Hilbert space of the theory, we use local quan- 
tities which transform covariantly under gauge transformations. Whereas the generator of 
gauge transformations in terms of Af is complicated by the noncovariant transformation 
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properties of the vector potential, when expressed in terms of gauge covariant variables it 
simply turns into a (color) rotation generator. Gauge covariant quantities can furthermore 
be contracted with themselves in color and lead to gauge invariant variables. In terms of 
such gauge invariant variables, Gauss' law, or gauge symmetry, becomes manifest. 

Nevertheless, not any choice of gauge covariant variables is appropriate. An appro- 
priate set of variables should describe the correct number of gauge invariant degrees of 
freedom at each point of space, and should also be free of ambiguities (such as, for in- 
stance, Wu-Yang ambiguities, where several gauge unrelated vector potentials may lead to 
the same color magnetic field 0,0). In I., the set of gauge covariant variables, uf, we 
have chosen to define is given by the following differential equations: 

e^Djut = e ijk {d jU l + f abc A)ul) = . (1.1) 

The linear operator e % i k Dj = (Si) lk Dj, where S is the single gluon "spin" operator, plays 
a central role in our formulation. When Af is a pure gauge, the eigenvalues of this operator 
are ±p and 0; in this case the zero modes are the "longitudinal" gluons. In general, the zero 
mode wavefunctions replace the vector potential as the dynamical coordinate. Further, the 
remaining spectrum and eigenfunctions of the operator enter in the process of obtaining 
an expression for the electric field conjugate to Af in terms of the variables uf. It is clear 
that the spectrum is gauge invariant and that the wavefunctions uf transform as vectors 
under gauge changes. In I. and here we actually only consider in detail the SU (2) theory, 
j-abc _ £ abc (although in I. the extension to SU(N) is also partially treated). Because most 
of the details of calculation needed for our purposes have been spelled out in I., here we 
will rather present a brief summary of previous results. 

It turns out that for the SU(2) theory in canonical formalism, there is a natural 
symmetry under coordinate reparametrizations, which is respected by all commutators 
and basic formulas except for the Hamiltonian itself. We have purposefully maintained 
this symmetry in defining our new variables, so that a natural geometric picture arises 
as a guiding principle in the formalism at no extra cost. Under this reparametrization 
symmetry, the vector potential transforms as a covariant vector, while both the electric 
field E a% = —i5/5Af and the magnetic field, 

B ai = £ *jk( djA a + 1 £ abc A b A c^ ? ^ ^ 
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transform as contravariant vector densities. The canonical commutators between coordi- 
nates Af and momenta E h i also transform covariantly, and the Gauss law generator, 



g a = DiE ai = d l E ai + e abc A\E cl , (1.3) 

transforms as a scalar density. It is only because of this reparametrization covariance that 
we introduce the above seemingly peculiar placement of space indices. The only failure in 
reparametrization covariance comes in the Hamiltonian: 

H = - J d 3 x ( ° + g 2 s E ai E ai ) . (1.4) 

The integrand above is not a geometric scalar with the correct density weight, and the 
contraction in space indices is made with a Kronecker 8ij rather than with a metric tensor, 
that is, the Hamiltonian is "committed" to a flat space. 

It is straightforward to check that the definition of uf in (|1 . 1|) is identical to the stan- 



dard geometric equation defining the spin connection in terms of the Christoffel connection 
(or vice- versa): 

dj ut+e abc A b jU t -r>» =0, (1.5) 

where uf is a 3-bein, e ahc A\ is a spin connection, and r* fc is the Christoffel connection of 
the metric gij = ufuj. Requiring that uf transform as a vector under both gauge and 
reparametrization transformations then gives us a gauge and reparametrization covariant 
definition, and the above simple geometric picture. The "metric" tensor gr^ = ufu°- neatly 
organizes the six local gauge invariant degrees of freedom of the problem into a symmetric 
3x3 matrix, and the next task would be to write the Hamiltonian in terms of these 
variables. One can then prove that any gauge invariant functional of A® can be written 
as a function of gij only, and that any functional of gij is gauge invariant (cf. I.). This 
implements gauge invariance exactly. Further, it is also easy to include other color variables 
into the formalism when uf is used as the independent variable. 

Before we proceed to write down the gauge invariant, geometric expressions for the 
quantities of interest, we shall first observe that there are zero mode problems associated 
with the calculation of the electric field E ai = —iS/SAf. It is easy to see that, under 



the transformation of variables (|1 . 1|) , the Jacobian matrix 5Af/Suj involves the operator 
e ij Dj, which may have more than one remaining zero mode when the potential is not a 
pure gauge (where that is of course uf itself). It then seems there will be an indefinition in 
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expressing the electric field through a chain rule in terms of derivatives 8/8uf. Providing a 
careful treatment of this problem is one of the main goals of this paper. To establish that 
such a redundancy can be handled, a detailed treatment will be presented of an example 
where there remain in the presence of a large set of vector potentials, an infinite set of 
such zero modes. We shall for that purpose propose an infinitesimal deformation of ( |1 . 1| ) 
which, as we shall verify explicitly, resolves the difficulties with zero modes. In the sections 
that follow, we shall first write down all the relevant geometrical formulas in the presence 
of the deformation. Then, as an example of the treatment where this care is required, we 
will specialize to those gauge field configurations for which the associated geometries are 
3-spheres, where we shall be able to give the explicit expression for the electric field. In the 
limit in which the deformation is eliminated, we shall see that rather than an indefinition 
in the electric energy, there will be a restriction on the possible wavefunctionals describing 
states of the theory for such configurations. Finally, we shall give the expression for the 
electric field of a system of infinitely heavy point color sources immersed in these spherical 
configurations, using the formalism for introducing sources also presented in I.. 



In introducing a deformation to (|1 . 1| ) , again we are careful to preserve both 



reparametrization and gauge covariance, and we must verify that it indeed removes any 
zero mode ambiguities. We choose 



ljk D jU a k = pe ljk e ahc u b j u% , (1.6) 



where p is a small parameter with dimensions of mass. It is possible to find Af explicitly 
as a functional of uf. With manipulations similar to those found in I., one finds 

AiM=ri + { * mkB **>£?-M« ) . (1.7) 

One can already glean from the above why this eliminates zero mode problems: vari- 
ations Su, 8 A must satisfy 

e ljk (8 ac dj + e abc {A) - 2pu b j ))8u c k = -e^ k e ahc u)8A c k . (1.8) 

To obtain Su? in terms of 8Af, the operator acting on 8u® must be inverted. With p = 
this operator would be e l ^ k Dj itself, whose zero modes are the very solutions to (|1.1|) , while 
for p 7^ 0, on the other hand, possible zero modes of the operator on the l.h.s. of ( |1.8[ ) are 
clearly not solutions to (|1.6|) . In fact, the claim we shall make is that the operator in ( |1.8| ) 
has no zero modes for small enough nonzero p and can always be inverted, leading to an 
unambiguous definition of the electric field in the M-variables. We now proceed to present 
the relevant formulas in geometric variables. 
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2. Gauge Invariant Geometric Variables 



With the definition (|1.6|) , and using the gauge Ricci identity analogously to what was 
done in I., one can determine that the magnetic field, when expressed in terms of geometric 
variables, is 

B ai = y/g (G ij +p 2 g ij ) u a 3 . (2.1) 

Here Gij is the Einstein tensor of the metric g^, Gij = Rij — and throughout indices 

are raised and lowered with the metric g^ and its inverse g % K By ^fg we mean det u, which 
can in principle take on both positive and negative values. The gauge Bianchi identity can 
be worked out to be 

AB" = V9 i^iG ij ) u° = , (2.2) 

so that it implies the geometric Bianchi identity and vice-versa. We note that this is a 
nontrivial consistency check of our geometric picture. 

We now turn to Gauss' law and the electric field. If we define, following I., the gauge 
invariant tensor e lJ through 

= - «i°e* , (2-3) 

we can then verify that the Gauss law generator in geometric variables becomes 

%Qa = Di (n?) = V^ M i(^ eii ) = VgtfiVie^) +pu i* £ijke jk } (2 .4) 

where Vj is the geometric covariant derivative with the torsion term introduced by p ^ 0, 
and u ia is the matrix inverse of uf (or, what amounts to the same thing, Uj with the space 
index raised by g lJ ). Alternatively, it follows directly from (|1.6| ) that if one makes a gauge 
variation of Af, 

5A1 = -A5w a , (2.5) 

then 

£ ijk D (2 P ) ^ _ £ abc 5w b u cj = q ? ^ 

where the operator D^ 2p ^ contains in place of A the "potential" A — 2pu. Since eD^ 2p ^ has 
no zero modes it follows that, as expected, a gauge variation of A is equivalent to a gauge 
variation of u transforming as a vector so Q a is represented in terms of u by 



To calculate the electric field in the u- variables, we begin from the definition (lO) of 

e ij : 

e = u aJ —— = u aj —r- ■ — r , (2.8) 

where the dot stands for a generalized contraction including an integration over space. 
Variations in uf can be further separated into variations of the six gauge invariant degrees 
of freedom g^j and of three gauge degrees of freedom in the following simple way: first we 
split the summation over the index b by inserting the unit color matrix in the form 5 bc = 
u b m u cm . The resulting quantity u crn -J-cCai\. then be written as the sum of its symmetric 



and antisymmetric pieces in m and k. Finally, it is easy to see that these correspond, 
respectively, to variations in the metric gij and gauge variations. The final expression we 
arrive at is: 

^"/^^(^SS^^C^^S^^) • (2 - 9) 

where i 

ga _ £ abc u b 

8/811% is again the Gauss law generator, now in terms of the u- variables, 
and we make explicit the space integration. 

We now need to write the Jacobian matrix 8u/SA in geometric form. To do so, we 
start by considering the following eigenvalue problem: 

£ iik( 5 ac dj + £ a bc{A b _ pu b ))wA e = ^/gX A W A ia , (2.10) 

where again indices are raised with the inverse metric g^ = u ia u ja '. We note that, by 
definition, one solution to the above with Xa = is uf itself. In the notation we are using, 
A is an index that labels all these eigenfunctions except the particular one given by uf. 
The operator above is real and symmetric, and we assume forms a complete 

orthonormal spectrum of real eigenfunctions for it. By orthonormality we mean 

J d s x ^u a lWA a 3 = 

J d 3 x y/ggftwAfwB* = 3V5 AB (2.11) 

J d 3 x y/gg ij uy> = 3 J d 3 x Jg = W , 

where V is the volume of the space described by g^ , and S A b is a Kronecker or Dirac delta 
depending on whether the spectrum is discrete or continuous. Because we will eventually 
concentrate on spherical geometries, we are only considering here configurations of finite 
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volume (that is the "dynamical" volume V, and not the volume of space, which is infinite). 
The generalization to infinite V should not entail further conceptual difficulty. 
If we now expand a generic variation 5u1 in terms of this complete set, 

5u?= V u? + J2vawa? , (2.12) 

A 

substitute this in (|1.8|) and dot it on the left with the same complete set ("dot" meaning 
an inner product with the measure *Jgg %:) ), we easily get the following relations: 

3Vpq = J d 3 x ^u ia 5A1 

r (2-13) 

wY^IabVb = - d 3 x e^e^utwA^At , 

B ^ 

where 

Iab = XaSab -^J d3x e^ k e ahc u a lWA h 3 w B l . (2.14) 
The origin of the zero mode problems alluded to above and their resolution through the 



p 7^ deformation of (LT) now become manifest: for p = 0, the first of eqs. (|2.13|) actually 
represents a constraint on the variations 5 A® for which a Suf can be found. This constraint 
is a direct consequence of the fact that ( |1.7| ) is homogeneous in uf for p = 0. For p ^ 0, this 
homogeneity is clearly broken, and there is no longer a constraint. Furthermore, essentially 
the same happens with the second set of equations in ( |2.13|) : for p = further constraints 
on SAf follow for each mode for which Xa = 0. Again, these are eliminated by taking 
p^O. 

From here, it is straightforward to write the Jacobian matrix in explicit form: 

Jj^) = zvp-^(y> Jb (yK(x) - w Y, I AB^ bcd <(y)^Bi(y)w A U x ) ■ (2.15) 

This can be expressed in an entirely geometric form by using the important fact that 
if WAi is a mode of ( |2.10|) , then the geometric modes za/ defined through 

waI = z Ai j u« (2.16) 

can be seen to be eigenmodes of the geometric curl operator £ ljfc Vj with the same eigen- 
values Xa- 

e ijk V j z Ak rn = ^gX A ZA irn ■ (2.17) 
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This leads to the fully geometric form we were seeking for the operator appearing in (gT 

^(x)^^it(y))=y/g(x) ^g ij (x)g km (y)+H j \Jx,y)-g i ^x)n s s krn (x,y) 

(2.18) 

where the Green's function Hij mn (x, y) is defined to be: 

Hijmn(x,y) = 7^^2zAij(x)lAB Z Bmn(y) ■ (2-19) 
AB 

Assembling all these results leads to the following electric geometric tensor acting on 
functionals 



e ij '(x)tf = / d 3 y 



' g ij (x)g mn (y) + {W l rnn (x,y)-^(x)H s s rnn (x,y)) 



3Vp 



Sgmn(y) 2 y/g(y) 



(2.20) 



From this expression, one may already observe that, independent of the geometry, there is 
always at least one divergence in the electric energy as p — > 0. We eliminate it by requiring 
that gauge invariant functionals &[u] be invariant under global rescalings of the metric, 
i.e., 

I d3v<iv) ^m =0 - (2 - 21) 

The Green's function Wy mn (x, y) may also have divergences in the limit p — > 0, which 
again have to be eliminated. Generally speaking, the higher the degree of symmetry of 
a certain geometry (which is determined by its Killing vectors), the larger the number of 
zero modes of the curl operator and, due to (|2.14|) , the larger the number of divergent 
terms in Tlij mn as p — > 0. In the following section we will work out and analyze the 
electric field for those compact geometries with the maximum number of Killing vectors, 
namely, spheres. Because they are maximally symmetric spaces, for spheres it is possible 
to find the spectrum of eV explicitly without too much difficulty, and therefore an explicit 
expression for Tiijmn as well. 

A final note on renormalization of divergences is in order here. It has been usual in 
the past literature on the subject to consider as the electric energy density expectation on 
a state \& the expression 

< *|(E-(x)) 2 |* >=- j [DA] , (2-22) 



with [DA] an appropriately denned integration measure. However, another way to define 
the electric energy expectation is (cf. I.) 



The former expression is inherently divergent due to the coincident points in the double 
functional derivative, and one must go to some lengths to properly define the operator. The 
latter expression, on the other hand, is easier to define and can be seen as an alternative 
prescription for the electric energy density. In our work we always use this second form. 
Of course, for systems with a finite number of degrees of freedom the two expressions are 
equivalent. Here, they differ formally by a total functional derivative. 

3. Spherical Configurations 

In order to get a clear picture of our restriction to spherical configurations, the first 
questions we will address are a) to what A® (x) configurations do spherical geometries 
correspond, and b) how much of the entire space of Af (x) do these geometries cover. 

The direct way of answering the first question is of course to take a configuration uf(x) 
describing a 3-sphere and substitute it in (|1.7|) to find Af(x). We will do so below for a 
particular metric on S3. There is, however, a more indirect but extremely economical way, 
based on the following reasoning: the geometry of a sphere is that of an Einstein space, for 
which Gij oc gij ; this implies, by ( |2.1|) , that the magnetic field for such geometries must be 
proportional to the matrix inverse of the 3-bein, B ai oc u ai . If we now consider the standard 
expression ( |1.2| ) for B ai as a function of Af, and the fact that pure gauge configurations, say 
Af, have vanishing magnetic field, it follows immediately that for configurations which are 
global scalings of a pure gauge, say Af(x) = kAf(x), k 7^ 1, the magnetic field will turn out 
to be proportional to B ai oc e % ^ k e ahc A h -A c k . But this is proportional to the matrix inverse of 
Af and thus, for spherical geometries, Af(x) must be proportional to uf(x). Closer scrutiny 
of this argument shows that it indeed holds and furnishes all the proportionality constants 
missing above. We now list a series of results that derive from the above reasoning. In 
what follows, we take A?(x) to be a pure gauge configuration, and a a real number 7^ 1. 




(2.23) 
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ii) Af(x) = 2(i-q) A-i(x) =>• uf(x) is an Einstein space, with Gij(x) = — ^gij(x) 

-V a 'i 9ij[X). 



6 ■ 

Af(x) — —uf(x) is pure gauge. 



Hi) uf(x) is an Einstein space, with Gij(x) = —^gij(x) = — p 2 ^ 1 a ^ 9ij( x ) ==> A® (x) 



iv) Af(x) = A^{x) — ^uf(x) is pure gauge =>- (A 9 )^(x) = ^(u 9 )°;(x) for some gauge 
transform (A 9 ,u g ) of (A,u). 

Thus, up to gauge transformations, the circle of implications above flows freely in both 
directions. Moreover, we can now also answer the second question as well: the space of all 
spherical geometries corresponds to the space of all vector potentials that are rescalings of 
all possible pure gauges. More concretely, all possible pure gauges in SU (2) are spanned by 
three real functions £ a (x), and rescalings are spanned by one real number. ^ This is indeed 
expected and consistent with the geometric picture, since all possible 3-sphere metrics are 
spanned by one real number (the inverse radius of the sphere, in units of p) and three 
real functions y l (x) (coordinate reparametrizations of a reference metric). This should be 
contrasted with six real, local functions, which parametrize the physical, gauge invariant, 
Hilbert space of the theory, so that, roughly speaking, 3-spheres span half the dimensions 
of this space. It is also possible to study the case of noncompact maximally symmetric 
spaces, i.e., 3-hyperboloids, although we will not consider these geometries here. The 
vector potential can easily be found by substituting the appropriate bein uf (x) in (|1.7|) . 
The procedure described above for spheres can be extended to hyperboloids by taking a 
complex, which would lead to complex potentials. Insofar as the symmetries are concerned, 
complex potentials do not spoil any of the reasoning above; however, in order to have real 
vector potentials in the end would require a complex gauge transformation in iv) above. 
Altogether we know this is possible since it is not difficult to obtain a real vector potential 
for this case. 

To give a concrete example, we can consider the projective metric on S3 used below 
(cf. (3.16) and below for the coordinate conventions). It is not difficult to find that in the 
limit of interest to us, p — > 0, the associated gauge field configuration is 

Am = -2 . (3.1) 

a z + \xr 



t Incidentally, the special case of a = 1 can be treated separately, and is seen to correspond to a 
flat geometry where, in an appropriate gauge, u1(x) — 6? and the vector potential is A2(x) = p8f . 
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Here, a -1 is the radius of the sphere. Two features of these configurations are worthwhile 
noting: first, if we take a 2 negative, we simply get the hyperbolic configurations men- 
tioned above; unlike spherical configurations, they have singularities at finite x. Secondly, 
these configurations correspond precisely to the spatial components of the instanton of 
Belavin et al, with a -1 the "size" of the instanton 0. Further, a closely related type 
of configuration for which such a geometrical picture is also readily available correspond 
to Wu-Yang monopoles || (which are in turn related to merons H). They are gotten by 
taking a -1 — > in the instanton configuration above, and multiplying it by a global factor 
|. The associated gauge invariant metric variable is gij = p 2 /\x\ 2 Sij, with p a constant 
parameter. This metric again describes a space of constant curvature, but closer inspection 
of its curvature invariants reveals that it actually corresponds to the space S2 x JR- It has 
been argued that (coherent states of) these magnetic monopoles are in fact one of the key 
ingredients underlying the color confinement mechanism in QCD ||. We will not pursue 
such a geometry further in this paper, but we wanted nonetheless to illustrate the point 
that our geometric setting fits very nicely with specific gauge field configurations that are 
deemed to be important for the dynamics of Yang-Mills theory. 

It is also curious to note that for spherical configurations, and again in the limit p — > 0, 
two of our fundamental equations are identical to the two equations defining the classical 
phase space of d = 3 gravity in the presence of a cosmological constant. For uf seen as a 
dreibein and uf 3 = e abc A^ seen as a spin connection, Eq. ( |1.1|) states that the connection 
is torsion free, and (|2.1| ) restricted to spheres states that the curvature built out of the 
spin connection is proportional to the inverse dreibein. These are just the equations of 
motion of d = 3 Einstein-Hilbert gravity with a cosmological constant (related to p 2 in 
( |2.1| )) or, equivalently, of an 50(4) Chern-Simons action with the gauge field being given 



by a combination of both and uf 3 [fLOf . This analogy, of course, does not go any further 



since it makes no mention of the electric field or of the specific form of the Hamiltonian. 

We now present our results on the spectrum of the curl operator ||. The most 
straightforward way we found of calculating this spectrum was by writing ansdtze for the 
modes ZAij based on the scalar, vector and 2-tensor eigenmodes of the Laplacian on the 
3-sphere and, after exhausting all possibilities for covariant ansdtze, verifying that the 
resulting modes ZAij satisfy a completeness relation. The eigenmodes we find are: 

1. Exact zero modes. 

z Aij = (Z N ) ij = (ViVj + ^g^yN ,N = 0,2,3,... . (3.2) 

a VKr -3) (oft - 1) 
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Here and in what follows, lvJj = N(N + 2) is the spectrum of the Laplacian acting on 
scalars on the 3-sphere, yN are the associated eigenmodes (the hyperspherical harmonics) 
normalized to 3V, and a is the inverse radius of the sphere (with scalar curvature R/6 = 
a 2 ), so that 

(V 2 + a 2 c4)yjv = 0. (3.3) 

The index A represents the three quantum numbers (N, £, m) labelling these modes. Be- 
cause I and m do not affect the spectrum, but only its degeneracy, they are of secondary 
importance here, and will be omitted and understood whenever possible. The normaliza- 
tion is such that geometrical modes satisfy the normalization conditions following from 
( [2. 11|) . For all these modes, Xa = 0. We also note that iV = 1 is missing: it vanishes 
identically due to the structure of the operator (VV + a 2 g). Moreover, the N = term 
will also be eliminated from the calculation of the Green's function in (|2.19| ) because it 
corresponds to the metric zero mode gij, which is to be treated separately (cf. discussion 
above (fmp). 

2. Scalar-based modes. These are nonzero modes based on the spectrum of the Laplacian 
acting on scalars: 

ZAij = (S Na )ij = — ~, = ( ViVj + a 2 u 2 N gij - aaJu% - l^V fc J y N . (3.4) 

2a 2 u N y/uf f - 1 V v V9 J 

Here, a = ±, and the associated eigenvalues are Xa = ^Na = cta^/ufj — 1, N = 1, 2, 3, — 

3. Vector-based modes. These modes are based on the spectrum of the Laplacian acting 
on vectors, 

(V 2 + a 2 (u 2 N -l))(v%) t = 0, (3.5) 
with V l (f^r)i = 0. Such vectors are also eigenvectors of the curl operator, 



: ' ik Vi(v%) k =aaJu 2 N -l(v%) l , (3.6) 



y/9 

with a = ±. The vector-based modes are then given by 



z Aij — (^iV^)ij 



(3.7) 



with eigenvalues 



(3.8) 
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where /3 = ±, A^f = af3/ (a^u 2 N — 3), and iV = 2, 3, 4, . . . (the N = 1 modes built in this 
way vanish identically, similarly to what happens in the exact zero mode case). 

4- Tensor-based modes. These modes are based on the spectrum of the Laplacian acting 
on tensors, 

(V 2 + a 2 (u 2 N -2))(T Nct ) lJ =0 , (3.9) 

with V*(2jVa)ij = 0, (TNa)l = an d (Tjva)ij symmetric. They are given by the (Tjsr a )ij 
themselves, with a = ±, N = 2, 3,4 . . . and eigenvalues \Na = ota^ u 2 N + 1 = aa(N + 1). 
Again, the putative N = 1 mode vanishes identically, and the spectrum starts from N = 2. 

It is now simple (but lengthy!) to calculate the matrix Iab, through ( |2.14| ) (with 
the gauge modes w referred to geometric modes z through ( [2.16| )). If we organize the 
matrix into five sectors, T (for tensor-based modes), V (for vector-based modes), S+ (for 
scalar-based modes with a = +), S— (for scalar-based modes with a = — ), and Z (for zero 
modes), a simple structure emerges: the matrix is diagonal in the T and V sectors, and 
the only non-diagonal couplings appear between the Z, S+ and S— sectors. Furthermore, 
its subblocks are diagonal in each and all of its nonvanishing sectors (e.g., TT, VV, ZS+, 
ZS-, ZZ, etc.). This allows for an explicit inversion, even though the matrix is infinite. 
The nonvanishing entries of Iab i n the different sectors are: 



I Ntfjap = (aay^N + 1 +P)^MNS a f3 (3.10) 
^N^Ma'P' = ~Y (^J u N + 1 +0\/ u} N-^j <W<W<W , (3.11) 

(where Smn includes orthogonality in £, m as well) and 

(2 ujN^/uj'ij - 3 uj n ^/u% - 3 \ 

u, Ny /u%=3 f(^-l) 3 / 2 -^ -1 f (3 . 12) 

U, N yfi$=l -1 -fK-1) 3 / 2 -^/ 

in the Z, S+, S— sectors, where the first, second and third row (or column) refers to, 
respectively, Z, S+ and S— sectors. Each entry represents an infinite diagonal matrix and 
because of this, inversion can be accomplished by simply inverting the 3x3 matrix. This 
inverse is 



a 2 



-2 

c4 



2(p 2 - a 2 )i 

(P + £,n)^n a /^at ~ 3 -2£ N - pu 2 N p(u 2 N - 2) 

\(p-£n)unt/u'n ~ 3 p(u}%-2) 2£ N -pv 2 N 



(3.13) 
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with £ 



N 



1. 



These results can finally be substituted in (|2.19|) in order to calculate / Hij mn (x,y). 
The TT and VV contributions can easily be gleaned from (|3 .10 ) and ( 3.11 ), while the 
contribution from the Z, S+, S— sectors is quite lengthy. In fact, a series of simplifications 
take place, and the final result is: 

1 



n 



ijmn 



P 



(ViVj- + a z giJ ) x iy m W n + a 2 g mn ) y G 3 (x, y) + 



2a 2 (p 2 -a 2 ) 
a 



2a 2 (p 2 — a 2 ) 



2a 2 p 



(ViVj 



+ 



y/9 

V9 



ae r 



V9 



(V m V n ) 



G (x,y) 



G (x,y) + ... . 

(3.14) 

The dots represent the TT and VV contributions, and Go and G3 are the Green's functions 
for the operators V 2 and (V 2 + 3a 2 ), respectively, acting on scalars on the sphere: 

n , n, _ 1 yN(x)y N {y) 

N ~ X (3.15) 
J_ yN(x)y N (y) 



G3(x : y) 



3V ^ a 2 (u 

N=2 v 



3) 



Since (V 2 + 3a 2 ) has zero modes given by yu m (x), it does not strictly speaking have a 
Green's function; what we mean by the above is of course the Green's function on the 
subspace of functions on £3 that is orthogonal to this zero mode. In fact, both Green's 
functions above are also lacking the trivial yo = const, mode in their spectral sum. This 
will be reflected in the differential equations they satisfy. 

It is in fact possible, with some effort, to find closed expressions for these propagators. 
We present them in what follows and briefly describe how they are gotten since these are 
useful in the calculations envisaged in Sec. 4. We use the standard projective metric on 



g i j(x 1 ,x 2 ,x 3 ) 



(3.16) 



[1 + a 2 \x\ 2 ) 2 13 ' 

where Xi, x 2 , X3 are projective coordinates, with range —00 to 00, \x\ 
a is the inverse radius of the sphere. From this, the Laplacian acting on scalars on S3 can 
be built, leading to the following differential equations for Go and G3: 

—o^,, s 5(x — y) 1 / \ / \ 8(x — y) 1 
V °^ y) = ~ -^gW + W^ yM = ~^gW + V 

n i (3.17) 



(V 2 +3a 2 )G 3 (a:,y) = - 



S(x - y) 



+ 



^7 yN£m(x)y N em{y) 



N=0 £=0 m= 
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where V = 27r 2 /a 3 is the volume of the sphere. The extra terms on the r.h.s. represent the 
lack of completeness of these Green's functions, as alluded to above. We will not go into 
the lengthy details of calculation, but rather just present the final result for Go and G3: 

G (x, y) = ±( ^Z? _ ^1=] sin -i ) « (3 . 18) 

G 3 (x, y) = ^ r \ Ji - 8 VV /1 - V 2 ) (1 - -(sin" 1 V - cos" 1 77)) + t^(<V + 1) 



16tt \ rj^i-^ /v ' / TT 24tt 3 



(3.19) 



where 



_ ad a\x-y\ 

V = T = /1 t 21 if (3 - 20) 

is one half the chordal distance d between the points x and y in units of a -1 . 

We are in fact interested in the p — > limit of ( ^.14j ) , and there are a few important 
features to note regarding this limit. Firstly, in the TT and VV sectors, the p —>■ limit 
is perfectly smooth; in particular, we would have obtained the same result had we taken 
p = from the beginning. Secondly, we find in the ZZ sector another 1/p divergence 
exactly like the one associated with global scalings (cf. (|2.2G| ) and ( |2.21|) ). This divergence 



appears in the first term in (|3.14|) , and it will likewise entail a constraint on finite energy 
physical wavefunctionals. What this constraint is can be seen from the term g l -''H s s mn in 
( |2.20|) : the trace in the first two indices leads to the operator (V 2 + 3a 2 ) acting on G3, and 
this leads to three types of terms, as can be seen from ( |3.17| ). The first term is a 5(x — y), 
and vanishing of this term leads to the constraint 



Iran } 



Sgmn(y) 



= (3.21) 

3ij=sphere 



in order not to have a divergence in the limit p — ► 0. The second term contains yo(x)yo(y) 
const .. and the second operator (VV + a 2 g) acting on it kills it because of the constraint 
fl2.21|) . The third type of term contains yu m (x)yiem(y) ■, and these vanish automatically 
under the action of the second operator (VV + a?g). Thus, only the first term leads to a 
constraint, ( |3.21j ), on physical wavefunctions.This constraint must be satisfied by physical 
wavefunctionals in order for their energy to be finite in the limit p — > O.This is again a 
result we would obtain directly in the p = case from a similar requirement of finiteness 



of the electric energy. Finally, we note that the last term in ( |3.14|) above is finite and 
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nonvanishing in the p — > limit. This result is different from what one would obtain by 
treating the p = case directly. The correct result is the one presented here, since it 
properly takes into account the mixing of the zero and scalar modes. 

In the p — > limit and with the above finiteness constraints in place, the Green's 
function 7i; jmn is 

Hijmn{x,y) = 4)] \ — / ! [(TN+)ij(x)(T N+ ) mn (y) - (T N -) tj (x)(T N -) mn (y)] 

SV N ( a V W N- 1 

2 

+ : , 2 =- — , 2 =: [(VN++)ij(x)(V N++ ) mn (y) - (V N - + ) ij (x)(V N - + ) mn (y)} 

a W W N + 1 + V U N ~ 3 ) 
+ , , 2 ==T / 2 ^ [(VN+-)ij(x)(V N+ -) mn (y) - (V N —)ij(x)(V N —) m n(y)] 



1 

+ 



(ViV . )B f £^^j ^ + X {VmVn)l 



2a 2 

while its trace reduces to a single term 



G (x,y) , 

(3.22) 



K- mn ( X ,y) = -^(^Y S -^. (3.23) 



2a 2 V v 7 ^ / Vd( x ) 

From now on we assume this limit unless stated otherwise. We now have an explicit 
expression for both electric and magnetic energy densities. Although the electric energy 
is still in a rather unwieldy form, it is possible already to make an important observation 
regarding the vacuum state of the theory: spherical configurations introduce a scale into 
the problem, as any other explicit configuration would. Such a scale must be dynamically 
determined and, although we will not perform such a calculation here, we can already 
observe that this scale a enters the magnetic energy with a positive power and the electric 
energy in negative powers. This will cause the ground state wave functional to fall rapidly 
for large amplitude magnetic densities which vary slowly in space. At the same time, it will 
also become small for low amplitude magnetic energy densities with slow spatial variation. 
The correspondingly reduced fluctuations in the magnetic energy density will presumably 
fill the role of what is meant by a magnetic "condensate" . If this is correct then one should 
be able to get at least a semi-quantitative estimate of the long range color electric fields 
produced by static sources in such surroundings. To do this we suggest looking at the 
electric field by dropping terms which should mainly be associated with short scales. 
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It would also be possible to perform manipulations in the vector and tensor sectors 
similarly to what has been done for the scalar and zero sectors, in order to simplify their 
respective contributions to Hij mn - However, as we shall argue, it is these latter sectors 
(i.e., the VVVGo term in (|3.22|) ) which will give rise to the main contribution to the 
potential between color sources, and therefore one may in fact drop the TT and VV terms 
in a first approximation. To identify particular terms in TCijmn that lead to large electric 
energy densities, we must look both for eigenvalues Xa such that 1/ A a becomes large, and 
for modes which do not oscillate much, since highly oscillating modes cannot contribute to 
long distance effects. The third term in (|3.22| ) (the last VV term) does have asymptotically 
large inverse eigenvalues as N — > oo; however, these are also associated to highly oscillatory 
modes. On the other hand, all theslowly oscillating T and V modes do not have inverse 
eigenvalues that become asymptotically large. The only term satisfying both conditions 
we are seeking is the last one, associated to the scalar and zero sectors, and therefore it 
is reasonable to keep only this term as a leading approximation. Naturally, our formalism 
automatically guarantees, as announced in Sec. 1, that this represents a gauge invariant 
approximation to the dynamics. 

4. Static Point Color Sources 

We now consider the energy density of infinitely heavy point color sources immersed 
in the Yang-Mills configurations associated to spherical geometries. The formalism for 
introducing point color sources has been developed in I.. 

Introducing color sources at isolated points in space entails a local modification of 
Gauss' law only at these points. Then, rather than introducing an additional set of vari- 
ables at every point in space, one can accomodate these isolated inhomogeneities in Gauss' 
law by simply considering wavefunctionals that carry the appropriate representation for 
each source, but that are still functionals of only. To be specific, let us consider, for 
instance, the insertion of two sources at points X\ and x^- The generalization to more 
sources is entirely trivial, but we consider here this specific case for clarity of presentation. 
Then, wavefunctionals describing states of this system should take the form 
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(4.1) 



where a is an index in some SU(2) representation transforming at point^i and f3 likewise, 
but transforming at x%. The modification in Gauss' law is 



Q a {x) -> G a {x) aa >pf3> = G a {x)5 aa >5pf3> + A. a aa ,5pi3>5{x - x±) + 6 aa >App,6(x - x 2 ) , (4.2) 

where A a are the appropriate SU(2) generators. Again, to be specific, we consider the two 
sources to be in the fundamental representation, in which case the A a are proportional to 
the Pauli matrices a a . The statement of gauge invariance becomes 

Q a {x) aaip ^ al p,[u1} =0 . (4.3) 

At this point it may not be clear whether or how one can build a color singlet wave- 
functional, satisfying the local constraint Ql.2p, exactly at the locations x\ and x 2 of the 
sources, since at each of these points the total color has contributions coming only from 
the combination of "integer spin" variables in the adjoint representation (the u$), and a 
"half-integer spin" variable (a or f3) coming from the source, which is in the fundamental 
representation. This turns out to be possible because there are sufficient variables uf in 
order to build a half-integer spin representation of SU(2) at x\ and at x 2 with these vari- 
ables alone, even though they are in the adjoint representation. The way this is done is by 
realizing that u", and u$ form three vectors, and thus comprise nine degrees of freedom 
at the point x\ or x 2 . While six of these degrees of freedom, u^Uj (i.e., the gauge invariant 
ones), give three lengths and three angles with which to uniquely define a tetrahedron, 
the three remaining degrees of freedom can be used to define the Euler angles uniquely 
fixing the orientation of the tetrahedron in color space. We then make use of the fact that 
it is possible to build half-integer spin representations of SU(2) with three Euler angles. 
With the appropriate transformation of variables from uf to Euler angles, the angular 
momentum operator becomes precisely Q a , and the appropriate eigenfunctions are the 
Wigner -functions, with j = \ for angular momentum one-half. Equations (|4.2[) , (|4.3|) 



then express the fact that under the usual addition of angular momentum in quantum 
mechanics, the wavefunctional at X\ and at x 2 is a singlet, of total angular momentum 
zero, built out of two spin one-half representations. Because this procedure is to be done 
at the isolated points X\ and x 2 , that is, because of the delta functions 5(x — x\) and 
5(x — £2), the wavefunctional, besides being a functional of u^(x) everywhere, must now 
also be a regular function of the variables w"(xi) and ii"(x2): 

* = vI> a/3 K;<(x 1 ),<(x 2 )] , (4.4) 
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so that the functional differentiation in Q a at those points becomes a regular derivative, 
and automatically incorporates the two delta functions. 

With the introduction of sources, the change in magnetic energy is due not to any 
modification in ( |2.1[ ) but rather to the constraints on \& engendered by ( [4. 3| ) . For the 
electric energy, on the other hand, there are modifications coming from the Gauss law 



term in (|2.20| ) , which would otherwise be absent for gauge invariant functionals. Neglecting 
the TT and VV contributions to r Kij rnn in fl3.22|) , the contribution of the sources to the 
geometric electric tensor is 



2 source 



d 3 y 



V9 



G (x,y) 



-mnp 



y/g(y) 



(v m v r 



Go(x,y) 



(4.5) 



which simplifies to 



2a 2 e™(x)v%] = [(V l V j rVt 1 G (x,x 1 )]A k (x 1 )^+g lJ (x)-^(V x k 1 S(x - x 1 ))A k (x 1 )^> 

(4.6) 

plus an identical contribution at X2- Here, Ak(x) = — iu%(x)A a , and we have ommited the 
577(2) indices. The action of Q a on W has been such as to satisfy (|4.3| ). 

One could now calculate an expression for the "potential" associated with static 
sources just as one evaluates the static Coulomb energy in the abelian gauge theory. It is 
clear that as well as being more complicated, the static potential is a function of the gauge 
field configuration and hence even in an approximation of the Born-Oppenheimer type it 
must be averaged with the ground state wavefunctional of the gauge field. Here we have 
obtained an explicit expression only for those gauge field configurations which are related 
to each other by GL(3) transformations of equal curvature geometries. We shall postpone 
a more complete discussion, and an explicit evaluation of the "potential" associated with 
the electric field in (|4.6| ) for a later publication. 



5. Conclusions 

In this paper we have pursued further the formalism developed in I., where a set of local 
gauge invariant variables were introduced to describe the physical Hilbert space of Yang- 
Mills theory. We have chosen to do this in a Hamiltonian, fixed-time formalism because 
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there one can identify the subset of the full gauge group that truly acts as a quantum 
mechanical symmetry of the theory, and one can implement it in a manifest and exact 
way. We have furthermore showed that the present treatment allows for approximations 
to the dynamics that do not spoil this exact gauge symmetry. 

As a first step towards concrete calculations, we have worked out in detail the expres- 
sion for the electric and magnetic energies for those gauge field configurations corresponding 
to spherical geometries. This has furnished indications of the mechanism through which 
a dynamically determined scale enters the theory and leads to a nonvanishing magnetic 
energy density of the vacuum and a mass gap. We have also indicated explicitly how such a 
geometry is related to instanton configurations, and how magnetic monopole configurations 
also correspond to a fairly simple, constant curvature space - S2 x M - in our geometric 
formulation of the theory. Moreover, for spherical configurations, we have studied particu- 
lar terms in the electric field energy in the presence of heavy point sources that lead to the 
main contribution to the potential for these sources. We have also identified the manner 
in which exact gauge symmetry is maintained locally in the presence of half-integral spin 
sources, whereby one must construct half-integer spin representations from the gauge field 
variables in order to construct total angular momentum zero from the addition to the color 
sources. 

Our calculations are by no means complete, and a number of important issues must 
still be considered: for instance, we have not studied the Jacobian determinant appearing 
in the measure after the change of variables, det \8A/8u\, and we have not considered the 
effects of renormalization. We also expect infrared effects to appear once noncompact 
geometrical configurations are considered, and these must be properly treated. Such issues 
would form an integral part of a more detailed computation of, for instance, the potential 
energy between two static color sources. A more detailed study of the S2 x JR geometry 
would also be of interest. All such computations are part of our plans for future work. 
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